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Abstract 

Exclusive production of and J/'0 mesons in e'’'p collisions has been studied with 
the ZEUS detector in the kinematic range 0.25 < < 50 GeV^, 20 < VE < 167 GeV 

for the data and 2 < < 40 GeV^, 50 < VE < 150 GeV for the J/'0 data. Cross 

sections for exclusive p^ and J/ijj production have been measured as a function of 
Q^, W and t. The spin-density matrix elements Tqq, r\_i and Re r^Q have been 
determined for exclusive p^ production as well as Tqq and for exclusive J/V' 
production. The results are discussed in the context of theoretical models invoking 
soft and hard phenomena. 
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1 Introduction 


We report measurements of exclusive electroproduction of and J/il) mesons at electron- 
proton centre-of-mass energy of a/s = 300 GeV using the ZEUS detector at HERA. The 
reaction ep —> eVp, where V stands for a vector meson (p°, 0, J/^), is a rich source of 
information on soft and hard diffractive processes as well as on the hadronic properties of 
the virtual photon [Q. 

Exclusive photoproduction of light vector mesons (p°, uj and 0) has been studied in a wide 
range of the photon-proton centre-of-mass energy W both in hxed target experiments p| 
and at HERA For fU > 10 GeV these reactions display features characteristic 

of a soft diffractive process: s-channel helicity conservation, cross sections rising weakly 
with W and a steep exponential t dependence, where t is the squared four-momentum 
transfer at the proton vertex. Such processes are well described within the framework 
of Regge phenomenology and the Vector-Meson Dominance model (VMD) [0], where 
exclusive vector-meson (VM) production at high energies is assumed to proceed via the 
exchange of a Pomeron trajectory as shown in Fig. |l]a. In this approach, the W and t 
dependences are coupled: 


du^'P 

d|f| 


oc e 


-bo\t\ 


W' 


2X 2{a{t)-l) 


( 1 ) 


where a{t) = q;(0)-|-q:T, while bo and Wq are process-dependent constants. Fits to hadron- 
hadron scattering data and photoproduction data give a(0) = 1.08 and a' = 0.25 GeV^ [Q. 
The slope of the t distribution depends on the energy as b = bo + 2a'/W q) (often 
referred to as “shrinkage”), while the effective power of a dependence of the cross 
section (after integrating over f) is 5 ~ 4(a(0) — l—a'/b). Typically, values of 6 ~ 10 GeV“^ 
are found in the photoproduction of light VMs, leading to 5 ~ 0.22, in agreement with 
measurements. This approach fails to describe the recently measured energy dependence 
of the cross section for elastic J/0 photoproduction at HERA The measured slope 

of the |t| distribution, bj/^ ~ 5 GeV“^, leads to a prediction of h ~ 0.14, in contrast to 
the measured value of S ~ 0.9. 

Exclusive VM electroproduction at high values of has been studied in fixed target 
experiments 000 and at HERA |T^, |T^. The measurements indicate that the rise 
of the cross section with W is stronger than that expected from Regge theory although 
there are large uncertainties in the experimental data|]^. The dependence of the cross 
section can be described by with 2 < n < 2.5 and the |t| dependence has a slope b 
between 4 and 8 GeV“^. The vector mesons are found to be produced predominantly in 
the helicity 0 state, whereas in photoproduction the production is mainly in the helicity 
±1 states. 

In models based on perturbative QGD (pQGD), 7 *p —Vp scattering is viewed as a se¬ 
quence of events separated in time in the proton rest frame [0, as depicted in Fig. [pD. 
The steps are: the photon fluctuates into a qq state; the qq pair scatters on the proton 
target; and, the scattered qq pair turns into a vector meson. The interaction of the qq 


1 







pair with the proton is mediated in leading order by the exchange of two gluons in a 
colour singlet state. In this framework, the cross section is proportional to the square 
of the gluon density in the proton. The scale /i^, at which Og and the gluon density are 
evaluated, can depend on the mass of the vector meson My, on and on t. For J/ip 
production the scale jj? = + |f|]/4 ||^ has been proposed. In photoproduc¬ 

tion at small |t|, the scale would therefore have the value /i^ = 2.4 GeV^. At this scale, 
the gluon density at small x rises as xg{x) oc x' 0, yielding a W dependence of 


the cross section ct'^p oc significantly steeper than expected from VMD and Regge 

phenomenology. The calculation has been extended |0 and compared to HERA data. 
It was found that this process is indeed sensitive to the form of the gluon density in 
the proton. Other pQCD calculations have been performed within the framework of the 
Colour Dipole Model (CDM) where diffraction is viewed as the elastic scattering of 
a colour dipole of dehnite size off the target proton. 

At large Q^, the cross section is predicted to be dominated by longitudinally polarised 
virtual photons scattering into vector mesons of helicity state 0 0 0,0,0. The cross 
section, calculated in leading Og In ^ In ^ approximation |0 for vector mesons composed 
of light flavours, is 


dal^ 


d|f| 


t=o 




1 + 


ITT 


2 din a: 


xgix.Q"" 


( 2 ) 


where A is a constant which depends on the VM wavefunction. Here we discuss some of 
the expectations for the cross sections: 


The cross section contains a 1/Q® factor. However, the dependences of as and 
the gluon density also need to be taken into account. The effective dependence 
using the CTEQ3L gluon density function [0] and the leading-order form for as 
is approximately dcr^ ^/d|f| oc l/Q^, with a weak x dependence. The calculation 
presented in |]0| has been redone in leading Q;sln(Q^/A^) approximation |0. In 
this work, among other improvements, the Fermi motion of the quarks in the vector 
meson has been considered. The net effect is to reduce the steepness of the 
dependence. Precise measurements could therefore yield information on the wave- 
functions of the vector mesons. 


In the pQCD calculations, the t and W dependences are not coupled, so that no 
shrinkage is expected. A lack of shrinkage, along with a steep W dependence, 
indicates that the reaction is predominantly driven by perturbative processes. In 
such processes, the transverse size of the qq pair is small, and the slope is determined 
by the proton size, resulting in a value for b near 5 GeV~^ 


H- 


The cross section presented in Eq. is for longitudinally polarised photons. The 
authors of ref. |0| expect that this is the dominant contribution to the cross section 
in DIS. It has been argued |0] that the region of validity of the pQGD calculations 
is signalled by the predominance of VM production in the helicity zero state. A 
recent pQGD calculation for electroproduction |0], based on the production of 


2 






















light qq pairs and parton-hadron duality, gives an estimate of the transverse photon 
contribution to the 7 *p ^ p°p cross section. 


The interaction should be flavour-independent at sufficiently high scales. From the 
quark charges of the vector mesons and assuming a flavour-independent production 
mechanism, the exclusive production cross sections should be in the proportions 
9 : 1 : 2 : 8 for : ut : (j) : J/iIj- This expectation is badly broken at low Q^. The 
pQCD predictions change the ratio somewhat due to wavefunction effects, such that 
the relative contribution of heavier vector mesons is enhanced at high Q^. 


In this paper, we investigate the dependence of and J/il) production on the variables 
IF, and t. The and J/ijj mesons are identihed via their decay to two oppositely 
charged particles. Invariant masses are reconstructed under the assumption of dipion hnal 
states for the and dimuon hnal states for the J/^j. The decay angular distributions 
are also measured, and the helicity matrix elements extracted, yielding a measurement of 
i? = (t 2 ^ as a function of and IF. The results are compared to expectations from 

Regge theory as well as from pQCD models. The data are presented for production 
in the ranges 0.25 < < o.85 GeV^ (BPC p°) and 3 < < 50 GeV^ (DIS p°). The 

production of J/'0 mesons is investigated in the range 2 < < 40 GeV^ (DIS J/'ip)- 

The data discussed in this paper correspond to integrated luminosities of 6.0 pb“^ (DIS 
and J/'ip) and 3.8 pb“^ (BPG p°) collected in 1995. 


2 Experiment 


The measurements were performed at the DESY ep collider HERA using the ZEUS de¬ 
tector. In 1995 HERA operated at a proton energy of 820 GeV and a positron energy of 
27.5 GeV. A detailed description of the ZEUS detector can be found elsewhere 
main components used in this analysis are described below. 


m. The 


The high-resolution uranium-scintillator calorimeter GAL consists of three parts: for¬ 
ward [| (FGAL), barrel (BGAL) and rear (RGAL) calorimeters. Each part is subdivided 
transversely into towers which are segmented longitudinally into one electromagnetic sec¬ 
tion (EMG) and one (RGAL) or two (FGAL, BGAL) hadronic sections (HAG). The energy 
resolution of the calorimeter, determined in a test beam, is <JeIE = 0.18/V^ for electrons 
and (TeIE = O.db/A/E for hadrons, where E is expressed in GeV. 


Gharged-particle tracks are reconstructed and their momenta determined using the central 
(GTD) and rear tracking detectors (RTD) [^. The GTD is a cylindrical drift chamber 


operated in a magnetic held of 1.43 T produced by a superconducting solenoid. It consists 
of 72 cylindrical layers, organised in 9 superlayers covering the polar angular region 15° < 


^Throughout this paper the standard ZEUS right-handed coordinate system is used: the Z-axis points 
in the direction of the proton beam momentum (referred to as the forward direction) and the horizontal 
X-axis points towards the centre of HERA. The nominal interaction point is at X = U = Z = 0. 
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6 < 164°. The RTD is a set of planar drift chambers located at the rear of the CTD, 
covering the polar angle region 162° < 9 < 170°. 


The positions of positrons scattered at small angles with respect to the beam direction 
are determined in the beam pipe calorimeter (BPC) and the small-angle rear tracking 
detector (SRTD). 


The BPC is an electromagnetic sampling calorimeter consisting of 2 modules, placed one 
on each side of the beam pipe, 294 cm away from the nominal ep interaction point in 
the rear region. Each module is equipped with 26 tungsten plates (roughly 24 radiation 
lengths), separated by layers of scintillator hngers (strips) each 8 mm wide. The strips 
alternate in the horizontal and vertical orientation, providing two-dimensional position 
information. The energy and position resolutions for electrons in the BPC, measured in 
a test beam, were found to be 17 %/{E in GeV) and ~ 1 mm, respectively. 


The SRTD is attached to the front face of the RCAL. It consists of two planes of scintillator 
strips, 1 cm wide and 0.5 cm thick, arranged in orthogonal orientations and read out via 
optical hbres and photomultiplier tubes. It covers the region of 68 x 68 cm^ in X and Y 
with the exclusion of a 10 x 20 cm^ hole at the centre for the beam pipe. The SRTD has 
a position resolution of 0.3 cm. 


The luminosity is determined from the rate of the Bethe-Heitler process e’''p —> e’'' 7 p, 
where the high-energy photon is measured with a lead-scintillator calorimeter (LUMI) 
located aX Z = —107 m in the HERA tunnel downstream of the interaction point in the 


positron flight direction [K 


3 Kinematics and cross sections 


We will use the following kinematic variables to describe exclusive VM production (see 
Fig. i: k, k', P, P', q, the four-momenta of the incident positron, scattered positron, 
incident proton, scattered proton and virtual photon, respectively; = —q^ = —{k — k'Y, 
the negative four-momentum squared of the virtual photon; = {q + P)"^, the squared 
invariant mass of the photon-proton system; y = {P-q)/{P-k), the fraction of the positron 
energy transferred to the photon in the proton rest frame; and, t = [P — P')"^, the squared 
four-momentum transfer at the proton vertex. 

The kinematic variables were reconstructed with the so-called “constrained” method, 
which involves the momenta of the decay particles measured in the CTD and the polar 
and azimuthal angles of the scattered positron in the BPC {p^ at low only) or in the 
main ZEUS calorimeter and the SRTD (all high-Q^ events). Neglecting the transverse 
momentum of the outgoing proton with respect to its incoming momentum, the energy 
of the scattered positron can be expressed as 

Ee' ~ [2Ee - (Ey - Py)]/{1 - COS^e'), (3) 
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where is the energy of the incident positron, and py are the energy and longitudinal 
momentum of the VM, and 9^/ is the polar angle of the scattered positron. The values of 
and t were calculated according to 

= 2Ee>Ee{l + COs6*e'), (4) 

\t\ = iPe' + Pvf + iPe' + Pv)\ (5) 

where , p^, and py, py are the X and Y components of the momentum of the scattered 
positron and VM. The variable y was calculated according to the expression 

y = {Ey - py)/2Ee (6) 

and Bjorken-x was evaluated using the relation = sxy, where s is the squared ep centre 
of mass energy. The kinematic ranges covered by the data are shown in Fig. |^. 


In the Born approximation, the positron-proton cross section can be expressed in terms 
of the transverse, and longitudinal, virtual photoproduction cross sections as 


dV"P 

dpdg2 


rT(y,Q^)(a^^ + eal p ). 


(7) 


where T^ is the transverse photon flux and e is the ratio of longitudinal and transverse 
photon fluxes, given by e = 2(1 —p)/(l-l-(l —y)^). In the kinematic range of this analysis, 
the value of e varies from 0.94 to 1.0. The transverse photon flux is p2 


g l + (l-y)^ 
2tt yQ^ 


( 8 ) 


where a is the hne-structure constant. The virtual photon-proton cross section cr'>' p = 
cr^ ^ can be used to evaluate the total exclusive cross section, a^ot = ctt ^ + o'l 


through the relation: 


,7"P _ 
Otot — 


1 +/? 
1 -I- ei? 


a 


7 P 


(9) 


where R = ctl is the ratio of the cross sections for longitudinal and transverse 

photons. 


At given values of W and the exclusive production and decay of VMs is described by 
three angles: <I>h - the angle between the VM production plane and the positron scattering 
plane; 0h and 0h ^ the polar and azimuthal angles of the positively charged decay particle 
in the s-channel helicity frame, in which the spin quantisation axis is dehned along the 
VM direction in the photon-proton centre-of-mass system. The angular distribution as a 
function of these three angles, IF(cos6*h, 0h, 4>h), is described by the spin-density ma¬ 
trix elements, pf^, where i, fc=-l,0,l and by convention a=0,l,2,4,5,6 for an unpolarised 
(or transversely polarised) electron beam 


The superscripts denote the decomposition 
of the spin-density matrix into contributions from the photon polarisation states: unpo¬ 
larised transverse photons (0), linearly polarised transverse photons (1,2), longitudinally 
polarised photons (4), and from the interference of longitudinal and transverse amplitudes 
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(5,6). For given values of W and Q^, the polarisation parameter e is constant, so that the 
contributions from and pf/. cannot be distinguished. The decay angular distribution 
can therefore be expressed in terms of linear combinations of the density matrix elements, 
and as 


r 


04 

ik 


r 


a 

ik 


p^k + ^Rptk 

1 + eR 
' Pfk 
1 + ei?’ 

'/R pfk 

. 1 + ei?’ 


a = 1,2 

a = 5, 6. 


( 10 ) 

( 11 ) 


Under the assumption of s-channel helicity conservation (SCHC), the angular distribution 
for the decay of the meson depends on only two angles, 6*h and '^h = 0h — "hh, and is 
characterised by three independent parameters, rQQ,r\_^ and Re rfg, as 


hF(cos0h,^h) = ^[^(1-’^oo) + - 1) cos^^h 

+ sin^ 6*h cos — 2Y^e(l + e)Re(r^Q) sin 26*h cos -^h] • (12) 


The ratio R 



can be determined from the polar angle distribution via 


R 


„04 

00 


e 1 — r 


04 • 
00 


(13) 


The additional assumption of natural-parity exchange in the t-channel reduces the number 
of independent parameters to two. The polar and azimuthal angle distributions are related 
via 


r 


1 

1-1 


1 

2 


(l-ro"o), 


(14) 


independently of the value of R. 
cesses at low energy . 


These relations were found to hold for diffractive pro- 


Statistical considerations limited our helicity analysis of the J/'tp sample to the one¬ 
dimensional distributions in 6'h and 0h, integrated over $h. For the decay to spin-1/2 
fermions, the above assumptions of SCHC and natural-parity exchange in the f-channel 
yield distributions sensitive to Tqq and according to 


lU(cos0h) 


3 

8 


1 + ^00 + (1 - 3roo) cos^ 0^ 


(15) 


W (cos 0h) = 1 + ’"l-l cos 20h 

271 L 


( 16 ) 


A value for R can be extracted from the polar angle distribution using Eq. 13. 
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4 Event selection 


The online event selection is done with a three-level trigger system The exclusive 
reaction ep —> ep°p at low (BPC) was selected at the first trigger level by the re¬ 
quirement of an energy deposit in the BPC of more than 6 GeV and at least one track 
candidate in the CTD. The DIS and J/'0 trigger at the hrst level performed an initial 
identification of the scattered electron in the main calorimeter by looking for isolated 
energy deposits. 


At the second trigger level general timing cuts were applied, along with a cut on the 
quantity E —pz = J2iiEi —pzi)-, where Ei = \Jp\^i + PY,i + Pz,i denotes the energy in the 
i-th calorimeter cell. The latter cut rejected background from photoproduction events. 
The BPC trigger included a restriction on the number of tracks in the CTD. 


At the third trigger level, requirements specihc to the exclusive reaction were imposed. 
These requirements were similar to those applied offline and included a vertex cut, a 
limit on the maximum number of tracks reconstructed by the CTD and a restriction 
on the maximum energy in the inner rings of the FCAL (rejecting events with proton 
dissociation). In the case of the DIS and J/'ijj events a positron candidate in the 
CAL was required and a fiducial cut on the positron position close to the rear beam pipe 
performed. 

In the offline selection of the exclusive and J/'ijj candidate events, the following further 
requirements were imposed: 


The energy of the scattered positron was required to be greater than 20 GeV if 
measured in the BPC and greater than 5 (DIS p°) or 8 GeV {J/'ijj) if measured in 
the uranium calorimeter. Positron identification in the latter two analyses used an 
algorithm based on a neural network |^. The efficiency was greater than 96%. 
In the case of the BPG, cuts were imposed on the deposited energy, shower width, 
timing and the BPG fiducial region. With these cuts, the probability that a selected 
particle is a positron exceeds 99% 


E — pz > 40 GeV. This cut, applied in both DIS analyses, excluded events requiring 
large radiative corrections. 

The Z coordinate of the interaction vertex was required to be within ±50 cm of the 
nominal interaction point. 


• In addition to the scattered positron the presence of two oppositely charged tracks 
was required, each associated with the reconstructed vertex, and each with pseudo- 
rapidityQ \r]\ less than 1.75 and transverse momentum greater than 150 MeV. These 
cuts excluded regions of low efficiency and poor momentum resolution in the tracking 
detectors. 

^The pseudorapidity rj is defined as ry = —ln[tan(|)]. 
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A match between each of the aforementioned tracks and an energy deposit in the 
calorimeter was reqnired. Energy deposits not associated with tracks or the positron 
were required to be less than than 300 MeV (elasticity cut), using a matching pro¬ 
cedure developed for this analysis . 


In addition, the following cuts were applied to select kinematic regions of high acceptance. 
The BPC analysis was limited to the region 0.25 < < 0.85 GeV^ and 20 < kh < 

90 GeV. The DIS analysis was restricted to the kinematic region 3 < < 50 GeV^ 

and 32 < kb < 167 GeV. For the cross section calculation, only events in the 7r’''7r“ mass 
interval 0.6 < < 1.2 GeV and with |f| < 0.6 GeV^ were taken in both analyses. 

For the J/il) analysis, cuts of 2 < < 40 GeV^ and 50 < kV < 150 GeV were applied. 

Only events within the mass interval 2 < Mi+i- < 4 GeV were accepted, where Mi+i- was 
calculated using the muon mass for the J/ijj decay products. 

The above selection procedure yielded 5462 events in the BPG sample, 3039 events in 
the DIS /9° sample and 213 events in the J/ijj sample. 


5 Acceptance corrections 


In the BPG analysis, a dedicated Monte Garlo generator based on the JETSET p9 
package was used to evaluate the acceptance and resolution of the ZEUS detector. The 
simulation of exclusive p° production was based on the VMD model and Regge phe¬ 
nomenology. Events were generated in the region 0.15 < < 1.1 GeV^, and 15 < kP < 

110 GeV. The effective W and t dependences of the cross section were parameterised 
as (T^ot oc l/(l-|-(5^/M^)^-’^®, a^ot oc kP°'^^ and d(T®P/d|f| = exp(—6|f| -1-ct^) (6 = 9 GeV“^, 
c = 2 GeV”'^), respectively. Decay angular distributions were generated assuming SGHG. 
A sample of events was generated using HERAGLES [40| in order to evaluate the magni¬ 
tude of radiative corrections in the BPG p° data. For the selected events they were found 
not to exceed 2% for any data point and to be consistent with zero within statistical 
errors. A 2% error was thus included in the normalisation uncertainty. 


In the DIS p° analysis, a dedicated program interfaced to HERAGLES was used 
to evaluate the acceptance and resolution associated with the “constrained” method of 
reconstruction. The cross section for exclusive p° production was parameterised in terms 
of (Tl ^ and over the entire kU and range covered by the data. At high initial- 
state radiation (ISR) introduces not only an overall correction to the cross section but also 
signihcantly distorts the distributions of certain kinematic variables. In the “constrained” 
method, ISR leads to migration of events along lines of constant W towards higher values 
of Q^. Moreover, it leads to additional and biased smearing of the reconstructed value 
of t. (For a cut on E — pz of 40 GeV, smearing due to ISR produces a decrease of the t 
slope by 5%). Final-state radiation does not introduce a signihcant error, as the radiated 
photon is usually well contained within the calorimeter cluster of the scattered positron. 


In the DIS J/^|J analysis the Monte Garlo program DIPSI [^, based on the model of 
Ryskin |^, was used. In this model, the exchanged photon fluctuates into a cc pair 




























which subsequently interacts with a gluon ladder emitted by the incident proton, and 
SCHC is assumed. The longitudinal and transverse cross sections are related by Ul ^ = 
where Mj/^ is the mass of vector meson J/^l]. For the gluon density in 


the proton the MRS A' parameterisation was used |Q. Events were generated assuming 
an exponential t distribution exp(—6|t|) with 6 = 5 GeV~^. The same method and tools as 
in the BPC analysis were used to calculate radiative corrections. Their magnitude was 
found not to exceed 4% for any data point and this value was included in the normalisation 
uncertainty. 

Two other generators, PYTHIA and EPSOFT, were used for determination of the back¬ 
ground from processes in which the proton dissociates. For the DIS and J/'ip analyses, 
the EPSOFT [^] Monte Carlo, developed in the framework of HERWIG [^, was used. 
It was assumed that the cross section for the reaction 7 *p —> VN, where N denotes the 
hadronic final state originating from the dissociated proton, is of the form 


dV^*P^VN ^ ida^T^Vp /^^pp^pN ^daPP- 
Mj.J 71 /r2 - ■' ' I - - - -o / 


^PP 


(17) 


where the ratio 


do-pp- 

dtdM, 


pN , 
N 


do-PP-*PP 


dt 


is obtained from fits to pp data [44 


The PYTHIA 

generator [|^ was used for the BPC analysis. A cross section of the form 

e“^'*'Esd(AfN)/AfN is assumed in PYTHIA, with 6 = 6o -|- 2a' In (bF^/M^), 6o = 2.8 GeV“^ 
and a' = 0.25 GeV^, corresponding to an effective 6 ~ 5 GeV~^ in the kinematic region of 
our results. The function Fsd{M^) enhances the cross section in the low mass resonance 
region and suppresses the production of very large masses |^. A fit to the generated 
spectrum for 10 < < 200 GeV^ with a function of the type 1/M^ gives n = 2.2. The 

effect of the functions Fsd{M-!<s) and 6 = 6(Mn) on the spectrum is thus consistent with 
the result n = 2.24 ± 0.03 measured for the diffractive dissociation of the proton in pp 
collisions ^ . 


In all three analyses the generated events were processed through the same chain of 
selection and reconstruction procedures as the data, accounting for trigger as well as 
detector efficiencies and smearing effects in the ZEUS detector. The distributions of 
generated variables were reweighted so as to reproduce the measured distributions after 
reconstruction. Gorrections for the data, evaluated on the basis of the Monte Garlo 
samples, were calculated independently in each bin of any given variable. 


A comparison of data and MG simulation is presented in Fig. (Q^, W and Ee>), Fig. ^ 
{t] and pt) and Fig. ^ (cos^h and pJh)- The J/tp sample is restricted to the mass range 
2.85 < M\+\- < 3.25 GeV in order to reduce the contribution from hadron pairs. The 
dominant remaining background originates from the Bethe-Heitler process ep —e UR p, 
and this contribution is represented by the shaded areas in the respective figures. This 
process, where the lepton pairs are either electrons or a muons, proceeds via the fusion of a 
photon radiated by the incoming electron and one radiated by the proton. Single particle 
distributions are very sensitive to the correct simulation of the W, Q^, t and decay angle 
variables. As an example, the transverse momentum distribution of the decay pion in 
the DIS sample (Fig. ^ displays a two-peak structure with maxima positioned around 
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~ 0.3 and ~ 1.8 GeV^ (this shape is less distinct in the case of the BPC p°s). Since at 
large valnes of the mesons are predominantly prodnced in the helicity zero state, 
one of the decay pions is emitted along the direction of the while the other one is 
approximately at rest in the 7 *p centre-of-mass frame. This conhguration resnlts in the 
Pt spectrnm, measured in the laboratory frame, shown in Fig. The measured and 
simulated spectra of cos6*h and V’h (for |cos6*h| < 0.5 and |cos6*h| > 0.5) are shown in 
Fig. 1^. The polar and azimuthal angular distributions are strongly correlated and the 
observed agreement between the measured and the simulated distributions was obtained 
via careful tuning of the simulation. 


6 Background 


After applying the selection criteria described earlier, the data still contain contributions 
from various background processes: 


• Proton-dissociative vector-meson production, ep —*• eVN, where iV is a state of mass 
Mn into which the proton diffractively dissociates. 

• Elastic production of uj and 0 mesons (for the analyses) and of ip' mesons (for 
the J/ip analysis). 

• Photon diffractive dissociation, ep —eXp and ep —eXN, in which the photon 
diffractively dissociates into a state X and the proton either remains intact or dis¬ 
sociates. 

• Bethe-Heitler production of e+e“ and pairs. 

• Beam-gas interactions. 


The main source of background consists of events with the proton diffractively dissociating 
into hadrons. Some of the hadrons deposit energy around the beam pipe in the FCAL 
but a fraction escape detection. The contribution by this process to the observed yields 
was estimated by using exclusive VM events with an energy deposit of at least 0.4 GeV 
in the FGAL. (Gontamination of this sample by DIS events is negligible.) 

Fig. I shows the ratios of the W, cos6'h and t distributions for FGAL-tagged events 
to those for all events in the BPG p° sample. The fraction of FGAL-tagged events is 
approximately independent of W and cos6*h. However, a signihcant dependence on 
t is observed. The latter is expected as a consequence of the different t dependences 
of the cross sections for elastic and proton-dissociative reactions. The same conclusions 
can be reached for the FGAL-tagged events from the DIS p° and J/ip samples. Under 
the assumption that a tag in the FGAL does not affect the shape of the acceptance as 
a function of W and cos^h (as indicated by PYTHIA and EPSOFT), this result 
suggests that proton-dissociative and elastic vector-meson production have the same 
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W and cos6*h distributions. This supports the hypothesis of factorisation of diffractive 
and in DIS [p0|] . 


vertices 


Similar conclusions were reached earlier for production at ~ 0 © 


In the BPC and DIS analyses the proton-dissociative background evaluation was per¬ 
formed by comparing the number of events with energy deposited in FCAL in the data 
and the proton dissociative EPSOFT or PYTHIA samples. Specifically, the number of 
residual proton-dissociative events in the data with FCAL energy smaller than the thresh¬ 
old Eq (1 GeV for the BPC analysis and 0.4 GeV for the DIS analysis) was estimated 
as 


y^DATA 

'^pdiss 


J ^ pdiss 

\ N{Efcal > Eq) 



x{iV(EFCAL>^0)}°^^^, 


where Apdiss is the fraction of elastic events passing the final cuts. A total of 160 (64) 
FCAL-tagged events were used for the BPC (DIS p°) analysis. In the BPC case, the 
additional requirement kF > 50 GeV was also imposed, in order to reduce the contribution 
by nondiffractive events. The overall contamination integrated up to |t| = 0.6 GeV^ was 
estimated to be (23 ± 8)% for the BPG sample and for the DIS sample, 

where the errors represent the total statistical and systematic uncertainties. In the DIS 
J/'ijj analysis, the contamination was found to be (2llg°)%, consistent with the values 
found in the BPG and DIS analyses. 


Gontamination from elastic production of uj and 0 mesons in the analyses was estimated 
by Monte Garlo studies to be less than 2%. Gontamination from tjj' production in the 
J/'ijj analysis is (4 ± 1)% |^. The Bethe-Heitler contribution is approximately 15%; its 
size was estimated from the LPAIR Monte Garlo simulation . 


The photon-dissociative background was studied with a sample of events generated using 
PYTHIA. The events which pass the selection criteria of the present analyses have a 
flat distribution in up to about ~ 1.4 GeV. If all events at = 1.4 GeV are 
ascribed to this process, a 3% upper limit on the contamination from photon dissociation 
is deduced. A similar result is obtained if an extra constant term is added to the Breit- 
Wigner function used for the mass fit. 


A contamination of 1.5% from beam-gas events was deduced from event samples derived 
from unpaired electron and proton bunches, to which all the selection criteria described 
earlier were applied. 


All subsequent results are shown after subtraction of the contributions from proton- 
dissociative and (for the J/'ip analysis) Bethe-Heitler events. The estimates of the other 
backgrounds were included in the systematic uncertainties. 


7 Systematic uncertainties 

The systematic uncertainties are dominated by the uncertainties in the acceptance, the 
proton-dissociative background and the number of p° or J/ij) signal events. Table 
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summarises the various contributions to the uncertainties in the integrated cross section 
for the three analyses. 

In the following sections, whenever a result for a given quantity was obtained from a £t, 
it should be understood that the corresponding systematic uncertainty was determined 
by repeating the £t for each systematic check. The differences between the values of the 
quantity thus found and its nominal value were added in quadrature. 

The trigger efficiency and its uncertainty were estimated, whenever possible, by using 
samples of events selected by independent triggers. The model dependence was investi¬ 
gated by comparing the acceptances obtained with various Monte Carlo generators, or 
with the same generator but with different input parameters. In particular, the sensitivity 
to the IV and dependences of the cross section in the generator were studied, as well as 
the sensitivity to R. Various electron hnders were used to estimate the uncertainty due to 
the electron identihcation in the DIS analyses. In the BPC case a signihcant contribution 
originates from the uncertainty in the BPC alignment; its effect increases with decreasing 
Q^. The sensitivity to the cuts mostly reflects the effect of the CTD-CAL-matching and 
track quality requirements. 

The contribution due to the extraction of the number of signal events, reflecting the 
sensitivity to the mass htting procedure, has been discussed previously 0 for the 
analyses. In the J/ijj analysis, this uncertainty is dominated by the sensitivity to the 
shape used for the subtraction of the nonresonant background. 

The uncertainties in the luminosity, trigger efficiency, photon flux determination, a;, 0 and 
ip' backgrounds, photon dissociation (for p^s), proton dissociation, beam-gas contamina¬ 
tion, and J/tp decay branching ratios are treated as overall normalisation uncertainties. In 
addition, for the BPC and J/ip analyses, the contributions from the signal extraction 
procedure and from radiative corrections were included in the normalisation uncertainty. 
The normalisation uncertainties are for DIS for BPC p°, and for the 

7/0, dominated by proton dissociation and, for the BPC p°, the signal extraction. 


8 Results 

8.1 Mass distributions 

Acceptance-corrected differential distributions dN/dM^+T^- for BPC and DIS samples 
are shown in Fig. |^. The vr+vr” mass spectra deviate from the shape of a relativistic 
p-wave Breit-Wigner function. This effect may be explained by the interference between 
nonresonant and resonant tt+tt" production amplitudes [^]. The differential distributions 
dV/dMn-+jr- were htted in the range 0.6 < < 1.2 GeV, in several intervals, 

using a parameterisation based on the Soding model , which accounts for the effect of 
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interference according to [49 


diV/dM^^ 


A 




Mi, - Mi + ,M,r, 


+ B 


( 18 ) 


where Mp and Fp are the nominal mass and width of the meson, respectively; B is 
the nonresonant amplitude assumed to be constant and real and A is a normalisation 
constant. The values of the meson mass and width obtained by fitting Eq. are 
768±3(stat.) MeV and 152±6(stat.) MeV for the BPC sample and 762±3(stat.) MeV 
and 146 ± 7(stat.) MeV for the DIS sample. The ratio B/A decreases with Q^, as 
shown in Fig. 


The uncorrected differential distribution dN/dMi+i- for the DIS J/'ifj sample, shown in 
Fig. was fitted in the range 2 < M\+\- < 4 GeV with the sum of a signal function 
and an exponentially falling background. The former is a convolution of a Gaussian 
resolution function with the J/ifj mass spectrum obtained using the DIPSI Monte Garlo 
generator [Q including bremsstrahlung. No positive muon or electron identification was 
performed; the muon mass was used in calculating M\+\-. The main contributions to the 
background are from oppositely charged hadrons and from the Bethe-Heitler process. The 
htted value of the J/'0 mass is 3.114 ± 0.006(stat.) GeV and the width of the Gaussian 
resolution function is 26 MeV. Integrating the htted function in the above M\+\- range 
yields a signal of 97 ± 12 J/'0 mesons. 


8.2 Total cross sections 

The total cross sections for exclusive and J/'0 electroproduction, ep —>• eVp, were 
determined using the expression 

N ■ A 

cr(ep ^ eVp) = (19) 

where N is the number of events in data, A the overall acceptance, C the integrated 
luminosity and A the correction for the proton dissociation background. For the we 
quote the integrated cross sections for |f| < 0.6 GeV^ and for the invariant mass range 
2mjr < AFn-jr < Mp + SFp, where m.,^ is the mass of a charged pion, is the invariant 
mass of the two pions, Mp is the nominal mass and Pp is the width of the resonance 
at the nominal mass. In the DIS p^ analysis, we correct to the Born level. 

Total cross sections for exclusive and J/tp production, ep —> eVp and 7 *p —Vp, are 
given in Tables ^ and ^ The cross sections in each and W interval are quoted at 
values close to the weighted averages in the bins. The 7 *p cross sections were obtained 
from the ep cross sections using formulae (0)”(iD- They are insensitive to the value of R 
since e ~1. 


13 















8.3 dependence 


Fig. H shows the cross section for the process 7 *p ^ p°p as a function of Q^. The 
low-Q^ data from this analysis have been htted with two VMD-motivated functions and 
the corresponding curves are shown in the upper plot. A £t to the function cr^ot^ oc 
[1 + i?(Q^)]/(l + yielded M^s = 0.66 ± 0.05(stat.)±0.10(syst.) GeV (here R = 

(Tl '’/o't ^ was taken as measured in this analysis). A £t to a tot oc 1/(1 + Q'^/M^Y gives 
n = 1.75±0.10(stat.)±0.29(syst.) for the entire sample (lFo=50 GeV). 

The 7 *p cross sections are shown for the DIS data for hxed IF as a function of 
in the four lower plots. The cross section measurement at = 27(13) GeV^ and W = 
80(120) GeV has been translated to IF = 70(110) GeV using the IF dependence measured 
in this analysis. The data are consistent with a simple power law behaviour for > 
5 GeV^. Fitting the points at > 5 GeV^ with the form yields n = 2.07 ± 

0.22, 2.51 ± 0.15, 2.15 ± 0.31, 2.29 ± 0.18 for IF = 50, 70, 90,110 GeV, where the errors 
are statistical only. The systematic uncertainties are typically 0.05. The dependence 
is consistent with being independent of IF, and averaging the four values yields n = 
2.32 ± 0.10(stat.). 


The data from this analysis are compared to previous HERA measurements in Fig. 

The cross sections are quoted at the IF values used by the HI collaboration |^. A 


comparison at hxed IF entails smaller translation uncertainties than a comparison at 
hxed as the IF dependence is much weaker than the dependence. A ht to the 95 
ZEUS DIS data is shown to guide the eye. The results from this analysis are in excellent 
agreement with the previous ZEUS results The HI data are systematically lower 

than the ZEUS measurements by approximately 30 to 40%. 


A commonly adopted form for the dependence of the J/ip cross section is aYt oc 
1/(1 + . This form was htted to the 7 *p —>■ J/ijj p cross section shown in 

Fig. 1^. The curve on the figure represents this function htted to the two measured data 
points (evaluated at lFo=90 GeV) yielding ?7,=1.58±0.22(stat.)±0.09(syst.). This result 
is consistent with the HI measurement of n = 1.9±0.3(stat.) [|T^] . 


8.4 W dependence 


The measured cross section for exclusive production as a function of IF for Q^=0.47, 3.5, 
7, 13 and 27 GeV^ is presented in Fig. |T^ The curves show the results of hts to the data 
using the function . The results of the hts are given in Table The IF dependence 
of the p° production cross section at low values of (BPG p°) rises slowly with W\ 
5=0.12±0.03(stat.)±0.08(syst.) for Q^=0.47 GeV^. This result is consistent with the 
value b = 0.16 ± 0.06(stat.)]';o;i5(syst.) measured in photoproduction |^. Averaging the 
data in the range 3 < < 20 GeV^ yields b = 0.42 ± 0.12(stat.©syst.), which indicates 

that the IF dependence increases with Q^. In Fig. the ZEUS data are compared to 
results from the NMG [T^, E665 |T^, and HI experiments. The NMG, E665 and 
HI data points have been moved to coincide with the values of the present analysis. 
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This was done according to the Q'^ dependence reported by each of the experiments. The 
HI points at Q^=13 GeV^ were obtained by translating the cross sections measured at 
Q^=10 GeV^ and Q^=20 GeV^ and taking a weighted average. The NMG measurements 
were moved from Q^=Q.9 to Q^=7.0 GeV^ and from Q^=11.9 to Q^=13 GeV^, using values 
of R from the model of Martin, Ryskin and Teubner to evaluate The 


E665 measurements were moved from Q^=0.61 to Q^=0A7 GeV^ and from Q^=5.69 to 
Q2=3.5 GeV^ 


The cross section for exclusive J/ijj photoproduction, measured at HERA and at low 
energies [^, shows a rapid rise with PE, approximately as The measured W 

dependence at higher values of is also consistent with this behaviour, as can be seen 
in Fig. 1^ In this hgure, the HI data points were scaled from Q^=16 to 13 GeV^. The 
curves, drawn to guide the eye, display a 1E° ® dependence. 


8.5 Ratio of J/?/; and cross sections 

The values of the ratio of J/ijj and (7*p) cross sections, measured at (5^=3.5 and 
13 GeV^, are given in Table The correlated errors, which include those associated with 
the proton dissociation background subtraction, with the uncertainties in the trigger effi¬ 
ciency and with the uncertainty in the determination of the luminosity, do not contribute 
to the uncertainty in the ratio. The ratio increases with as can be seen in Fig. 


8.6 Differential cross sections dcr®P/d|t| 


The differential cross sections for exclusive production, d(T®P/d|t|, were measured in 
several and W intervals; they are shown in Fig. |^for the full data samples. The distri¬ 
butions were fitted with an exponential function of the form exp(—6|f|) for |f| < 0.3 GeV^. 
Since a linear exponent is not sufficient to describe the data at higher |t|, the quadratic 
form exp(—6|t| -|- cR) was also htted to both the BPG and DIS data for |d < 0.6 GeV^. 
The results from the linear hts are b = 8.5 ± 0.2(stat.)±0.5(syst.)±0.5(pdiss.) GeV“^ 


for the BPG sample (0.25 < < 0.85 GeV^ and 20 < IE < 90 GeV) and b = 

8.1 ± 0.6(stat.)±0.7(syst.)±0.7(pdiss.) GeV“^ for the DIS sample, which covers the 
kinematic region depicted in Fig. Since a major contribution to the systematic un¬ 
certainty arises from the uncertainty associated with subtracting the proton dissociation 
background, the error from this source is explicitly quoted. In order to illustrate the 
signihcance of the quadratic term in the exponent of the fitted function, the uncertainty 
due to the systematic error in the parameter c is indicated by a shaded band in Fig. IK 
Detailed results of fits in PE and intervals are summarised in Tables ^ and |^. 


The slope parameter 6 as a function of PE and is displayed in Fig. The results for 
BPG are consistent with a slow rise with PE. The DIS results are consistent with 
no PE dependence. Both results show signihcantly shallower slopes than that measured 
in photoproduction. The GDM calculation is shown for comparison. Its prediction of 
a decrease with is in reasonable agreement with the data. 
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The linear-exponent fit for the entire DIS J/'tp sample (2 < < 40 GeV^ and 55 < hh < 

125 GeV) in the range \t\ < 1 GeV^ yielded b = 5.1 ± l.l(stat.)±0.7(syst.) GeV“^, a result 
consistent with the value of 4.6 ± 0.4(stat.)lo!6(syst-) GeV“^ obtained in exclusive J/ip 


photoproduction 


8.7 Shrinkage of the diffractive peak 


Shrinkage of the diffractive peak was studied by reweighting iteratively the energy and b 
dependence in the Monte Garlo simulation according to 
da^^P 


d|t| 


oc exp 


+ c„f) ■ 


( 20 ) 


where Wq is a constant, t and W are the generated variables, and bo, a(0) and a' are 
the parameters tuned to the best agreement between the simulated and measured dis¬ 
tributions. The fit for the BPG data yielded Q!(0)=1.055±0.016(stat.)±0.019(syst.) and 
Q;'=0.19±0.09(stat.)±0.09(syst.) GeV“^, showing evidence for shrinkage, in agreement 
with theoretical predictions [^. A similar analysis performed using the DIS data gave 
an inconclusive result. 


8.8 Decay angular distributions 


The spin-density matrix elements, Tqq, r 


1 

1-1 


and Re rfg, were determined by a two- 
dimensional maximum-likelihood fit of Eq. to the cos 6^11 and i/lh distributions. The 
results are presented in Table The corresponding values of the ratio R = as a. 

function of Q/^ are displayed in the upper plot of Fig. The ratio was evaluated according 
to Eq. O. The results indicate that the ratio increases with the dependence steeper 
at lower Q^. At high values of the longitudinal cross section dominates. The solid 
line represents the result of a fit to the BPG data of the form R = kQ^, which yielded 
K = 0.81 ± 0.05(stat.)±0.06(syst.). The dashed line represents the results of the QGD- 
based calculation of ref. , which describes the data well. The lower plots of Fig. show 
R as a function of W for = 0.45 GeV^ and 6.2 GeV^. For = 0.45 GeV^, the data 
indicate a slow decrease of R with W. However, they are consistent with no dependence 
within two standard deviations. For = 6.2 GeV^, the measurements indicate a slow 
rise of R with W. The prediction of the model of ref. ||^ is in good agreement with these 
results. 


The values of the spin-density matrix elements Tqq and r\_^ satisfy Eq. |T^ within ex¬ 
perimental uncertainties, and are thus consistent with natural-parity exchange in the t 
channel. 


The J/ip spin-density matrix elements Tqq and were determined by one-dimensional 
fits of formulae (P!5[ ) and ( 0 )- The results for the entire kinematic region covered by the 
data, for which < Q‘^ >= 5.9 GeV^ and < W >= 97 GeV, are 

Tqq = 0.29 ± 0.19 (stat.) Igjg (syst.), r^ti = —0.04 ± 0.20 (stat.) ^o ;22 (syst.). 
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Using < e >= 0.99, a value of R of 0.411 q; 52 (statistical and systematic uncertainties 
added in quadrature) was extracted, significantly less than the values measured for the 
at similar Q^. 


8.9 Forward longitudinal cross sections 


In order to compare our results to pQCD calculations, we extract the forward longitudinal 
cross section according to 


daZ'’ 


d|f| 


t=o 


R 


1 + Rj 


,7’P 

^tot 5 


( 21 ) 


where |f|max is the upper limit on |f| for which the cross section was calculated. A com¬ 
parison of the measured and the predicted x dependence of the longitudinal production 
cross section at various values of is shown in Figs. |l^pl| . The shaded areas indicate 
normalisation uncertainties due to the proton dissociation background subtraction, the 
measured values of R and of the slope parameter b. (For the highest value, extrapo¬ 
lations of the R and b values were used.) 


In the model of Frankfurt, Koepf and Strikman the hard diffractive production of 
vector mesons by longitudinal photons is calculated in the leading-order approximation 
(cKs In ^) using leading-order parton distributions. Rescaling effects are accounted for by 
introducing an effective scale, at which the gluon density is evaluated. The curves 
in Fig. |l^, which use the ZEUS 94 next-to-leading-order (NLO) gluon density parameter- 
isation [^, show the degree to which the rescaling attenuates the dependence. Also 
shown are the effects of two assumptions concerning the wave function which result in 
different Fermi motion suppression factors as calculated by the authors. The assumption 
of a hard Fermi suppression attenuates the dependence, as does the rescaling, but with 
a different x dependence. The two effects are of comparable magnitude in the kinematic 
region covered by the data. The measurements indicate that the assumption of a hard 
Fermi suppression together with the Q^-rescaling results in an overcorrection. Clearly 
a quantitative understanding of the higher-order QCD corrections is necessary before 
information on the gluon density and on the p meson wave function can be extracted. 


The model of Martin, Ryskin and Teubner is based on the parton-hadron duality 
hypothesis, applied to the production of qq pairs. A comparison of the predictions using 
various gluon density functions to the measured x dependence of the forward longitudinal 
cross section is shown in Fig. The MRSA' MRSR2 [^, and ZEUS 94 NLO 
gluon density parameterisations lead to similar predictions, whereas the prediction using 
the GRV94 parameterisation is considerably higher. In the context of this model, the 
data are sufficiently precise to distinguish between GRV94 and the other parton density 
functions. 

Fig. ^ compares the calculations of the two models described above with that of Nemchik 
et al. 0. which is based on colour dipole BFKL phenomenology. Here, the model of 
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Frankfurt, Koepf, and Strikman uses the ZEUS 94 NLO gluon density function with 
rescaling and no hard Fermi suppression. The curves for the model of Martin, Ryskin, 
and Teubner represent a calculation which also employs the ZEUS 94 NLO gluon density 
parameterisation. With these choices, the models describe the data reasonably well, taking 
into account the normalisation uncertainties. Note that the normalisation uncertainty due 
to the uncertainty in b is largely independent for the various values and dominates at 
high Q^. The model of Nemchik et ah underestimates the cross section over the entire 
kinematic range investigated. 


9 Summary and conclusions 


We have studied the exclusive electroproduction of ep —>• ep°p, and J/V’ mesons, 
ep —> eJ/V' p, in the kinematic range 0.25 < < 50 GeV^, 20 < fU < 167 GeV for the 

data and 2 < < 40 GeV^, 50 < fU < 150 GeV for the J/'ip data. The results can be 

summarised as follows. 


• The 7r’''7r“ mass spectrum for exclusively produced mesons shows a deviation 
from the relativistic p-wave Breit-Wigner shape. This can be explained in terms 
of the interference between resonant and nonresonant production amplitudes. The 
relative contribution of the nonresonant amplitude is found to decrease with and 
becomes consistent with zero at ~20 GeV^. 


• The dependence of the 7 *p —>• cross section at low (0.25 < < 

0.85 GeV^) can be described by the function oc 1/(1 + Q'^/MpY with n = 
1.75 ±0.10 (stat.)±0.29(syst.). At higher values of the dependence can be htted 
with the function oc with the average htted value of n = 2.3 ± 0.1(stat.), 
essentially independent of W. For the DIS J/?/ sample the data are described by 
the function oc 1/(1 ± with n=1.58±0.22(stat.)±0.09(syst.). 


The W dependence of the 7 *p —>■ p^p cross section exhibits a slow rise with W 
at low values of Q^. Parameterising the cross section as a'lot oc yields the ht 
result (5=0.12±0.03(stat.)±0.08(syst.) for Q^=0.47 GeV^ (BPG p°). This value is 
consistent with that measured in photoproduction as well as with predictions based 
on soft pomeron exchange |^. The slope becomes steeper with increasing Q^. For 
3.5 < < 13 GeV^ the average value is 0.42 ± 0.12. This is less steep than the 


value of (5=0.92±0.14(stat.)±0.10(syst.) measured in J/xf) photoproduction [[TO 
The cross section for J/ip electroproduction has a W dependence consistent with 
the steep dependence found in photoproduction. 


• The ratio a{J/"ip)/a{p^) increases with but does not reach the flavour-symmetric 
expectation of 8/9 at = 13 GeV^. 

• The t distributions for exclusive production are well described by an exponential 
dependence dcr®P/d|t| oc for \t\ < 0.3 GeV^ with 6 ~ 8 GeV“^. The slope 
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decreases at larger values of |f|. The Colour Dipole Model [Q gives a reasonable 
description of the data. A lower value, 5.1 ± l.l(stat.)±0.7(syst.) GeV“^, has been 
obtained in exclusive J/ijj production in the kinematic region 2 < < 40 GeV^ for 

|t| < 1 GeV^. This result is compatible with that for J/ip photoproduction. 

• The measurements in the range 0.25 < < 0.85 GeV^ exhibit a W de¬ 

pendence in the \t\ distribution. This may be interpreted as due to the shrink¬ 
age of the diffractive peak, predicted in Regge theory. In this context, we hnd 
a(0)=1.055±0.016(stat.)±0.019(syst.) and a'=0.19±0.09(stat.)±0.09(syst.) GeV“^. 
Tests for shrinkage in the DIS p° sample were inconclusive. 


The ratio of the cross sections for longitudinal and transverse photons, R = 
increases with and shows a weak W dependence. For > 3 GeV^ these depen¬ 


dences are well reproduced by the model of Martin, Ryskin and Teubner [^]. For 
J/ip electroproduction, R ~ 0.4 at = 6 GeV^, 
for the p° meson. 


in contrast to the value of i? > 2 


• The measurements of the forward longitudinal cross section, dcr^ ^/d|t||t=o; for 
production have been compared to the results of calculations based on several pQGD 
models. The present level of accuracy in the measurements allows quantitative 
distinctions between the various calculations. 


In conclusion, our results for exclusive p° production show the range 0.25 < < 

50 GeV^ to be a transition region, where, as increases, the relative contribution of 
continuum production decreases, and the longitudinal contribution to the total 

cross section increases and becomes dominant. These trends encourage efforts to describe 
this process using the methods of perturbative QGD. 

Exclusive J/ip electroproduction is consistent with expectations from pQGD. The expo¬ 
nential slope of the |t| dependence is approximately 5 GeV“^ and the W dependence of 
the cross section is consistent with the steep rise observed in J/ip photoproduction. These 
dependences differ from those measured in electroproduction at ~ where the 

\t\ dependence is steeper and the W dependence shallower. We also hnd contrasting val¬ 
ues for i? in J/ip and electroproduction. Thus and My are shown to play dissimilar 
roles in setting the scale of the process. 
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Contribution from 

BPC 

DIS 

DIS J/ip 

Luminosity 

1.1* 

1.1* 

1.1* 

Acceptance: trigger efficiency 

5.5* 

<1* 

<1* 

Acceptance: model dependence 

1-4 

5* 

1 

Acceptance: electron identihcation 

3-10 

<1 

<1 

Acceptance: dependence on cuts 

2-10 

6 

+6 

-11 

Acceptance: photon flux determination 

1* 

1 

1* 

Procedure to extract the signal events 

10* 

1* 

5* 

Proton diss. background subtraction 

10* 

+7 = 1 = 
-12 

+ 11 = 1 = 

-13 

Elastic uj and 0 production 

1.6* 

<1* 

- 

Elastic xfj' production 

- 

- 

1* 

Photon diffractive dissociation 

+ 0 = 1 = 

-3 

+ 0 = 1 = 

-3 

- 

Radiative corrections 

2* 

1* 

4 * 

Beam-gas interactions 

+0 :)= 
- 1.5 

+0 =(= 
- 1.5 

+0 =(= 

- 1.5 

Branching ratio 

- 

- 

2.2* 


Table 1: Typical values of relative contributions (%) to the systematic uncertainty in 
the integrated cross sections presented in Tables ^ and The starred values are the 
contributions to the overall normalisation uncertainty in each of the three analyses. 
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w 

[GeV] 

mmm 

#evts 

lEo 

[GeV] 

— 

IBiBI 

CJ®P 

[nb] 

Cr7-P 

[ph] 

BPC 

20-90 

0.25-0.29 

0.29-0.33 

0.33-0.38 

0.38-0.45 

0.45-0.55 

0.55-0.85 

1074 

941 

857 

869 

784 

937 

51.1 

0.27 

0.31 

0.35 

0.41 

0.50 

0.69 

4.99±0.15±0.47 

3.98±0.14±0.47 

3.88±0.13±0.33 

4.34±0.15±0.56 

4.55±0.16±0.49 

7.28±0.23±1.06 

5.07±0.15±0.48 

4.64±0.16±0.55 

4.14±0.14±0.36 

3.86±0.13±0.49 

3.41±0.12±0.37 

2.51±0.08±0.37 

20-27 

27-35 

35-45 

45-55 

55-70 

70-90 

0.25-0.85 

955 

1024 

994 

897 

1018 

574 

23.4 

30.9 

39.9 

49.9 

62.4 

79.8 

0.47 

5.52±0.18±0.64 

4.75±0.15±0.59 

4.62±0.24±0.61 

4.28±0.15±0.56 

4.89±0.15±0.50 

4.72±0.20±0.56 

3.16±0.10±0.37 

3.16±0.10±0.40 

3.19±0.16±0.43 

3.74±0.13±0.48 

3.61±0.n±0.38 

3.44±0.14±0.40 

DIS 

32-40 

40-60 

60-80 

80-100 

3-5 

254 

492 

401 

318 

36 

50 

70 

90 

3.5 

0.141 iO.Olltaoio 

0.256 ± o.onl^ij}^ 

0.211 ± 0.016l[J;[j}4 
0.186 ±0.016l[j;[j?^ 

0.310 ± 0.025^^;^^^ 
0.318 ±0.02lt°;°J^ 
0.376 ± 0.027;[|;^25 
0.443 ± 0.036;[|;^^^ 

40-60 

60-80 

80-100 

100-120 

5-10 

380 

331 

234 

193 

50 

70 

90 

no 

7 

0.101 ±0.007+^;^^]] 
0.089 ± 0.008l[J;[j[]4 
0.066 ± 0.007+[J;[j[j^ 
0.058 ± 0.006+^ 

0.075 ± 0.005l^;^^| 
0.095 ± 0.009t[|;[][|^ 
0.094 ± 0.010;[|;[j?^ 
0.109 ± 0.012t°;°°2 

41-60 

60-80 

80-100 

100-140 

10-20 

106 

88 

72 

no 

50 

70 

90 

120 

13 

0.023 ± 0.002l^;^^^ 
0.019 ± 0.0021°;°°^ 
0.014 ± 0.002l[];[|^^ 
0.025 ± 0.003l^;[]^^ 

0.021 ± 0.002l^;^^^ 

0.024 ± 0.003;[|;[j[|^ 
0.025 ± 0.004;[|;^[|^ 
0.030 ± 0.004t°;°°J 

55-96 

96-125 

125-167 

20-50 

27 

17 

16 

80 

no 

150 

27 

0.006 ± 0.001+^;^^{ 
0.004 ± 0.00ll^;[|^2 
0.004 ± 0.00ll^;[j[jj 

0.0033 ± 0.00071^;^^^! 
0.0045 ± 0.0012+[[;[5[|}^ 
0.0053 ± 0.0015+[[;[5[|}^ 


Table 2: Exclusive p° production cross sections for |t| < 0.6 GeV^ in various and 
W intervals. The BPC p° cross sections are calculated for the invariant mass range 
2m.^ < < Mp + STp. The cross sections are given at Ql and ITo values assuming 

the and W dependence from this analysis. The uncertainties do not include the 
normalisation uncertainties, which are ti 4 % for the DIS p° sample and, for the BPC 
p° sample. 
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W [GeV] 

g" [GeV^] 

# events 

ITo [GeV] 

g^ [GeV^ 

[pb] 

[nb] 

50-100 

2-7 

31±7 

70 

3.5 

79±18li^ 

21±5±3 

100-150 

2-7 

20±7 

120 

3.5 

56T1911? 

29 ± 10 ± 6 

50-150 

2-7 


90 

3.5 

135±26l^^ 

25 ± 5 ± 4 

50-100 

7-40 

24±6 

70 

13.0 

30±7+| 

6±2± 1 

100-150 

7-40 

29±7 

120 

13.0 

39±9l^ 

17±4±3 

50-150 

7-40 


90 

13.0 

69±12l?o 

10±2±2 


Table 3: Exclusive J/ijj production cross sections in various and W intervals. Values 
are quoted at Ql and ITo. The first error is statistical, the second systematic. The sys¬ 
tematic uncertainties include the normalisation uncertainty of 'tUx added in quadrature. 


g^ [GeV^j 

5 

0.47 

0.12 ±0.03 ±0.08 

3.5 

0.40 ±0.12 ±0.12 

7.0 

0.45 ±0.15 ±0.07 

13.0 

0.41 ±0.19 ±0.10 

27.0 

0.76 ±0.55 ±0.60 


Table 4: The values of the parameter S obtained by htting the W dependence of 
for exclusive production with a function a^ot oc . The hrst error is statistical, the 
second systematic. 


g^ [GeV^j 

ITo [GeV] 


3.5 

13.0 

90 

90 

0.06 ±0.01 ±0.01 
0.43 ± O.lOlg;]]^ 


Table 5: The ratio oi J/'ip and cross sections measured at two values. The hrst 
error is statistical, the second systematic. 
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fPo [GeV] 

Ql |GeV=| 

[GeV] 

b [GeV-2] 

BPC pO 

47 

0.45 

0.6 < < 1.2 

8.5 ±0.2 ±0.5 ±0.5 

25 



7.6 ±0.6 ±0.8 ±0.5 

35 



8.8 ± 0.8 ± 1.2 ±0.5 

50 



8.8 ±0.5 ±0.4 ±0.4 

74 



9.0 ±0.6 ±0.7 ±0.7 

47 

0.33 


8.6 ±0.4 ±0.6 ±0.5 

47 

0.62 


8.3 ±0.6 ±0.9 ±0.5 

47 

0.45 

0.6 < < 0.7 

9.7 ±0.7 ±0.8 ±0.5 

47 


0.7 < < 0.8 

8.5 ±0.5 ±0.6 ±0.5 

47 


0.8 < < 1.2 

7.8 ±0.6 ±0.6 ±0.5 

DIS 

67 

6.2 

0.6 < < 1.2 

o 1 +0.6+0.3+0.7 
O--‘--0.6-0.7-0.4 

50 

6.2 


o 0+0.9 +0.4+0.7 
O-^-0.9 -0.4-0.4 

70 

6.2 


o /i+l.l +0.3+0.7 
0-^-1.1 -0.8-0.4 

90 

6.2 


’-7 /j+l.l +0.2+0.6 

' *^-1.0 -0.9-0.3 

67 

3.8 


7 4+0.8 +0.1+0.7 
'•^-0.8 -1.0-0.3 

67 

6.8 


o c+1.0 +0.4+0.7 
O-O-0.9 -0.6-0.4 

67 

13 


Q 't+2.0 +0.5+0.6 
' -1.8 -1.3-0.3 

102 

28 


/I /I+3.5 +3.7+0.5 
^•^-2.8 -1.2-0.3 

67 

6.2 

0.6 < Mtttt < 0.8 

Q q+0.8+0.2 +0.8 

^•'^-0.7-0.5 -0.4 

67 

6.2 

0.8 < < 1.2 

7 7+1.0+0.7 +0.6 
-0.9-1.1 -0.3 


Table 6: The values of the slope parameter b obtained by htting d(T®P/d|t| oc in 
the range |t| < 0.3 GeV^ in various Q^, W, and M^+^- ranges of the BPC and DIS 
samples. The hrst line of each of the BPC and DIS sections indicates the results of the 
£t to the full sample. The first error is statistical, the second systematic and the third is 
the uncertainty resulting from the subtraction of the proton dissociation background. 
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Wo [GeV] 

Ql IGeV^l 

b [GeV-2] 

c [GeV-^j 

47 

0.45 

9.5 ±0.3 ±0.6 ±0.5 

4.0 ±0.7 ±0.8 ±0.4 

67 

6.2 

9.5 ±0.8 ±1.1 ±0.9 

6.1±1.3±1.7±0.5 


Table 7: The values of the parameters b and c obtained by htting (T®P/d|t| oc 
in the range |t| < 0.6 GeV^ for the full BPC and DIS samples. The mass range 
0.6 < < 1.2 GeV was used. The first error is statistical, the second systematic 

and the third is the uncertainty resulting from the subtraction of the proton dissociation 
background. 




ITo [GeV] 

Q'^o [GeV^] 

'00 

27.5 

0.45 

0.30 ±0.02 ±0.03 

45 

0.45 

0.26 ±0.02 ±0.02 

71 

0.45 

0.23 ±0.02 ±0.02 

47 

0.33 

0.23 ±0.01 ±0.02 

47 

0.62 

0.32 ±0.01 ±0.02 

46 

6.2 

0.71 ±0.02 ±0.03 

67 

6.2 

0.68 ±0.02 ±0.04 

92 

6.2 

0.75 ±0.03 ±0.03 

67 

3.8 

0.68 ±0.02 ±0.02 

67 

6.8 

0.74 ±0.02 ±0.04 

67 

14.1 

0.76 ±0.03 ±0.06 

Wo [GeV] 

Q'i [GeV"] 

^i-i 

27.5 

0.45 

0.32 ±0.01 ±0.03 

45 

0.45 

0.35 ±0.01 ±0.02 

71 

0.45 

0.36 ±0.01 ±0.02 

47 

0.33 

0.36 ±0.01 ±0.02 

47 

0.62 

0.31 ±0.01 ±0.03 

46 

6.2 

0.12 ±0.02 ±0.02 

67 

6.2 

0.15 ±0.02 ±0.06 

92 

6.2 

0.10 ±0.02 ±0.04 

67 

3.8 

0.15 ±0.02 ±0.02 

67 

6.8 

0.11 ±0.02 ±0.02 

67 

14.1 

0.08 ±0.03 ±0.07 

Wo [GeV] 

<35 [GeV^ 

Re rfg 

27.5 

0.45 

0.150 ±0.005 ±0.010 

45 

0.45 

0.136 ±0.004 ±0.010 

71 

0.45 

0.132 ±0.005 ±0.010 

47 

0.33 

0.133 ±0.004 ±0.010 

47 

0.62 

0.156 ±0.004 ±0.010 

46 

6.2 

0.10 ±0.02 ±0.03 

67 

6.2 

0.11 ±0.02 ±0.03 

92 

6.2 

0.11 ±0.02 ±0.03 

67 

3.8 

0.11 ±0.02 ±0.04 

67 

6.8 

0.10 ±0.02 ±0.03 

67 

14.1 

0.09 ±0.03 ±0.03 


1-1 


and Re rfg determined nsing Eq. p!2 


Table 8: The spin-density matrix elements rQQ, r 
for varions valnes of W and (BPC and DIS samples). The first error is statistical, 
the second systematic. 
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Figure 1: Diagrams illustrating the exclusive electroproduction of vector mesons a) via 
pomeron exchange, and b) via exchange of a gluon pair. 
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ZEUS 95 




Figure 2: The kinematic regions covered by the and J/ijj data samples used for this 
analysis. Line types: solid - constant W; dashed - constant polar scattering angle of the 
vector meson (6*v); and, dotted - constant positron scattering angle {9^'). 
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ZEUS 95 


• DATA 
— MC 






0 20 40 

Q^[GeV^] 


100 200 
W[GeV] 


20 25 

Ee- [GeV] 


Figure 3: Comparison of the measured and Monte-Carlo-simulated distributions for 
Q^, W and Eg', the energy of the scattered positron. The shaded area indicates the 
contribution from the Bethe-Heitler process. 
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ZEUS 95 


• DATA 
— MC 



Figure 4: Comparison of the measured and Monte-Carlo-simulated distributions for 
pseudorapidity, rj, and the transverse momentum in the laboratory frame, pt, for the 
positively charged decay particle. The shaded area indicates the contribution from the 
Bethe-Heitler process. 
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ZEUS 95 


DATA 

MC 



400 


200 
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400 
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111 



r 







0 1 0 2.5 5 0 2.5 5 

cos(0h) for |cos(0,^)|<O.5 for |cos(0^)|>O.5 


Figure 5: Comparison of the measured and Monte-Carlo-simulated distributions of cos 
and V'h- The iph distributions are shown for two ranges of cos^h- The strong correlation 
of these variables is evident. The shaded area indicates the contribution from the Bethe- 
Heitler process. 
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ZEUS 95 




> 

a> 

O 

< 

m 



Q^[GeV^] 


Figure 7: Top: The acceptance-corrected differential mass distributions, normalised 
to unit area, for the BPC and DIS samples. Line types: solid - fit based on the 
Soding model ^^(cf. Eq. |^); dashed - contribution from the p-wave Breit-Wigner term; 
dotted - interference term; dash-dotted - background contribution. Bottom: ratio B/A 
(cf. Eq. |T^. The open point associated with an arrow indicates the value measured in 
photoproduction . 
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ZEUS 95 



M,+[GeV] 


Figure 8: The differential mass distribution for the DIS J/il) sample (not corrected 
for acceptance). The distribution includes contributions from both e’''e“ and pairs. 

Solid line - signal function (see text), dashed line - background; shaded area - contribution 
from the Bethe-Heitler process. 
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Figure 9: The cross sections, a^o^, for exclusive p° production as a function of for 
various values of W. Top figure: the curves represent hts to the low-Q^ (BPC) data using 
the functions (l + i?)/(l + Q^/Mgfj)^ (dotted line) and 1/(1+ Q^/M^)"' (dashed line). The 
open point with the horizontal arrow indicates the value measured in photoproduction 
Four bottom hgures: the solid lines represent a £t of the form oc for 

5 GeV^. Only statistical uncertainties are shown. 
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ZEUS 95 



6 10 20 30 6 10 20 30 


Q^[GeV2] Q2[GeV^] 

Figure 10: Comparison of the dependence of measured by ZEUS and HI. 
The ZEUS points were moved to the W values quoted for the HI measurements. The 
lines represent the results of hts to ZEUS 95 data. The error bars on the ZEUS points 
represent the quadratic sum of statistical and systematic uncertainties. The normalisation 
uncertainties in the ZEUS measurements are indicated by the shaded areas. 
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ZEUS 95 



0 2.5 5 7.5 10 12.5 15 17.5 20 

[GeV^] 


Figure 11: The cross section, for exclusive J/ip production as a function of Q^. 

Line types: solid - function a^ot cc 1/(1 + fitted to the data points at Q^=‘i 

and 13 GeV^; dashed - extrapolation of the fit result to Q^=0. Only statistical errors are 
shown. 
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(Y*p ^ pp) [nb] 


I 



Figure 12: The cross section, a^ot ■, exclusive production as a function of W for 
various values of Q^. Only statistical errors are shown. The lines represent the htted 
parameterisation oc W^. 
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* Low energy photoproduction ❖ E665 □ H1 93/94 
o ZEUS 94 photoproduction A NMC • ZEUS 95 


'a. 

Cl 

T 

Q. 



10 


10 


W [GeV] 


Figure 13: Comparison of cross sections, cr^ot^, for exclusive production, as a function 
of W for various values of Q^. The error bars represent statistical and systematic errors 
added in quadrature. The solid lines represent the £t results shown in Fig. |T^. The 
dashed line is the prediction by Donnachie and Landshoff |^. The overall normalisation 
uncertainties are shown as shaded bands for the NMC and ZEUS data points. The 


NMC 
(see text). 


E665 and HI [O data points were interpolated to the indicated values 
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o ZEUS 94 

• ZEUS 95 
□ HI 94 

* E401 
0 E516 

_J_^^_ I _^^^^_I_I_^_ 

10 10 ^ 

W [GeV] 


Figure 14: The cross section, a'lot ^ for exclusive J/V’ production as a function of W for 
various values of Q^. The error bars represent statistical and systematic uncertainties 
added in quadrature, including the normalisation uncertainty in the ZEUS measurements 
of -i 5 %- Measurements from the hxed target experiments, E401 and E516 are 
included for comparison. The lines, drawn to guide the eye, correspond to the cross 
section parameterisation oc 
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ZEUS 95 



Figure 15: The ratio of the 7 *p cross sections for exclusive J/'0 and production as a 
function of IT^. The dashed line indicates the flavour-symmetric expectation of 


8/9. The inner error bars represent statistical uncertainties; the outer error bars indicate 
the quadratic sum of statistical and systematic uncertainties. 
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ZEUS 95 



|t| [GeV^] 



|t| [GeV^] 


Figure 16; The differential cross section d(T®P/d|t| for (0.6 < < 1.2 GeV) and 

J/ijj production. Only statistical errors are shown. Line types: solid - £t of the function 
d(T®P/d|t| oc for |t| < 0.3 GeV^ (|t| < 1 GeV^ for dashed - d(T®P/d|t| oc 
for |t| < 0.6 GeV^; the shaded band indicates the systematic error resulting from the 
uncertainty in the parameter c. 
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Figure 17: The slope parameter 6 as a function of W and for exclusive 7r’''7r“ produc¬ 
tion (BPC and DIS samples) in the range 0.6 < <1.2 GeV. The open point with the 

horizontal arrow indicates the value measured in photoproduction |49]. The inner error 


bars represent statistical uncertainties; the outer error bars indicate the quadratic sum 
of statistical and systematic uncertainties. The shaded areas indicate additional normal¬ 
isation uncertainties due to the proton dissociation background subtraction. The dashed 
lines represent predictions of the Colour Dipole Model (CDM) of Nemchik et ah 
the corresponding and W values. 
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Figure 18 

^7*P 


R = aria^ 


The ratio of the cross sections for longitudinal and transverse photons, 
T for exclusive 7r’''7r“ production in the range 0.6 < < 1.2 GeV as a 

function of and IF, evaluated assuming SCHC. The inner error bars represent sta¬ 
tistical uncertainties; the outer error bars indicate the quadratic sum of statistical and 
systematic uncertainties. The solid line represents the result of a £t to the BPC data of 
the form R = which yielded k = 0.81 ± 0.05(stat.)±0.06(syst.). The dashed line is 
a prediction of the model by Martin, Ryskin and Teubner using the ZEUS 94 NLO 


parameterisation of the gluon density |56 
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Figure 19: The measured forward longitudinal cross section, da^ '’/d|t| ^ , as a function 

of X for the DIS sample. The inner error bars represent statistical uncertainties; the 
outer error bars indicate the quadratic sum of statistical and systematic uncertainties. The 
shaded areas indicate additional normalisation uncertainties due to the proton dissociation 
background subtraction as well as the measured values of the R = ratio and the 


slope parameter b. The curves show the predictions by Frankfurt, Koepf and Strikman 25 


using the ZEUS 94 NLO gluon parameterisation [^. The full and dashed lines show the 
result of the calculation assuming hard Fermi suppression with rescaling < Q^) 

and without rescaling = Q^). The dashed-dotted and dotted lines show the result 
assuming no hard Fermi suppression with and without rescaling. 
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Figure 20: The measured forward longitudinal cross section, da^ '’/d|t| ^ as a function 
of X for the DIS sample. The inner error bars represent statistical uncertainties; the 
outer error bars indicate the quadratic sum of statistical and systematic uncertainties. The 
shaded areas indicate additional normalisation uncertainties due to the proton dissociation 
background subtraction as well as the measured values of the R = ratio and the 

f-slope parameter b. The curves show the predictions by Martin, Ryskin and Teubner 
and correspond to various gluon parameterisations, indicated as follows: full lines - ZEUS 
94 NLO [Q, dashed lines - MRSA' |^, dashed-dotted lines - MRSR2 |5^, and dotted 
lines - GRV94 ^ 
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Figure 21: The measured forward longitudinal cross section, da'l ^/d\t\ ^ as a function 
of X for the DIS sample. The inner error bars represent statistical uncertainties; the 
outer error bars indicate the quadratic sum of statistical and systematic uncertainties. The 
shaded areas indicate additional normalisation uncertainties due to the proton dissociation 
background subtraction as well as the measured values of the R = ratio and the 

slope parameter b. The solid line shows the calculation by Martin, Ryskin and Teubner 
using the ZEUS 94 NLO gluon density parameterisation. The dashed line shows the result 
of a calculation using CDM [OT . The dashed-dotted line shows the prediction of Frankfurt, 


Koepf and Strikman ||^ using the ZEUS 94 NLO gluon parameterisation with rescaling 
(QeS < Q^) no hard Fermi suppression. 
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